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Based on a mesoscopic theory we investigate the non-equilibrium dynamics of a sheared nematic 
liquid, with the control parameter being the shear stress <r xy (rather than the usual shear rate, 
7). To this end we supplement the equations of motion for the orientational order parameters by 
an equation for 7, which then becomes time-dependent. Shearing the system from an isotropic 
state, the stress-controlled flow properties turn out to be essentially identical to those at fixed 7. 
Pronounced differences occur when the equilibrium state is nematic. Here, shearing at controlled 
7 yields several non-equilibrium transitions between different dynamic states, including chaotic 
regimes. The corresponding stress-controlled system has only one transition from a regular periodic 
into a stationary (shear-aligned) state. The position of this transition in the cr xy -7 plane turns out 
to be tunable by the delay time entering our control scheme for a xy . Moreover, a sudden change of 
the control method can stabilize the chaotic states appearing at fixed 7. 



I. INTRODUCTION 

Shearing a fluid of orientable, e.g. rod-like particles 
is a fundamental example of a soft-matter system driven 
far away from equilibrium. The complex dynamics of 
such systems, examples being worm-like or cylindrical 
micelles [JlH, (polymeric) liquid crystals, blood cells, or 
colloidal rod (e.g. tobacco virus) suspensions, has be- 
come a focus of many experimental investigations Q, 
computer simulations 0-0], and theoretical studies over 
the last years (for a recent review, see [§|). Within the 
isotropic high-temperature (or low-concentration) state, 
application of external shear typically yields either con- 
tinuous flow-alignment or a discontinuous parancmatic- 
nematic transition 0] , dependin g on how close the under- 
lying equilibrium transition is |10L ITlj . Corresponding 
rheological properties such as the "constitutive curve" 
formed by the shear stress function of the shear 

rate 7 are already non-linear, implying a non-Newtonian 
behavior of the viscosity. More spectacularly, shearing in- 
side the nematic (i.e., orientationally ordered, yet trans- 
lationally disordered) phase can induce various bifurca- 
tions, that is, non-equilibrium transitions between differ- 
ent dynamic " states" characterized by a specific periodic 
or non-periodic motion of the nematic director 0, [T2I — 
Il5j . Examples are the tumbling or wagging behavior ob- 
served both by theory (l2j and in experiments of sheared 
tobacco viruses @. Further issues receiving much at- 
tention are the appearance of rhcochaos (chaotic stress- 
strain curves and/or orientational dynamics) [l2l - fl4l fl~6l— 
[20j . and, sometimes combined with that, spontaneous 
spatial-symmetry breaking associated to shear-banding 
[U, [2l| . In the latter situation, the system separates into 
domains characterized by different local shear rates. 

Many of these fascinating phenomena can be described 
by mesoscopic theories focussing on the motion of some 
suitable, coarse-grained dynamic variables. An exam- 
ple is the approach developed by Hess [H, [23| and Doi 
pil |25| . which is based on an equation of motion for 
the tensorial orientational order parameter, that is, the 
alignment tensor a(i). A constitutive relation then yields 



the link from a(i) to the stress tensor er and, thus, the 
rheology of the system. The theory is highly non-linear 
due to the use of a fourth-order Landau-de Gennes free 
energy functional describing the relaxation of a towards 
equilibrium. As a consequence of this non-linearity, the 
(five-dimensional) dynamical system involving the inde- 
pendent components of a generates complex orientational 
dynamics and associated rheological properties, even if 
the analysis is restricted to simple (Couette) shear ge- 
ometries and to spatially homogeneous systems (for ex- 
tensions of the Hess-Doi approach towards inhomogc- 
neous systems see, e.g., [2a - l29j ). Another class of meso- 
scopic theories focusses directl y o n the shear stress <7 xy 
as a dynamic variable [1, Il9l. T20I [30l [3l| . an example 
being the non-local Johnson- Segalman model [IH, |32| |. 
These models are capable of describing, on a quite general 
level, complex rheological behavior such as shear-banding 
Q , a drawback being that they reflect the (orientational 
and/or translational) dynamics within the underlying liq- 
uid only indirectly. 

In the present study we are interested in both, rheo- 
logical properties and orientational motion, of a complex 
fluid composed of rod-like particles. So far, most the- 
oretical studies investigating the orientational dynamics 
take the shear rate 7 as an external driving parameter. 
Experiments, on the other hand, are often conducted 
at controlled shear stress (see, e.g., (HHHI); moreover, 
there are various devices where either stress or strain 
can be controlled. Thus, it is an interesting question to 
which extent these macroscopic conditions influence the 
observed dynamics. 

In the present study we explore the implications of con- 
trolled stress (as compared to controlled shear rate) via 
an extension of the mesoscopic equations for a(i) [12, [H| . 
In order to fix the (time-averaged) shear stress at an ex- 
ternally imposed value, we supplement the five equations 
of motion for the nematic order parameters by a further 
differential equation for 7. The approach of o~ xy towards 
its target value is governed by a control time, r g . This 
control scheme is inspired by the method used in exper- 
imental devices: the shear rate is adjusted such that the 
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desired value of the shear stress is approached. Some pre- 
liminary results for stress-controlled systems have been 
published in Ref. (36l |. In the present paper we under- 
take a more systematic study focussing on spatially ho- 
mogeneous states with only one (constant) non- vanishing 
component of the shear rate (and stress-) tensor. This 
restriction excludes, e. g., the appearance of shear bands 
fiUli; shear- induced layering [33] , and other effects 38 1 
such as the emergence of secondary flow. However, we 
consider an analysis of the homogeneous system as com- 
pulsory before the full hydrodynamics with spatially in- 
homogeneous shear rates is performed. Nonetheless, the 
present control approach could, in principle, be applied 
to inhomogeneous systems (and cross-couplings between 
deformation rate and stresses) without further ado. In- 
deed, the mcsoscopic theory on which our study is based, 
is fully capable of treating inhomogeneous effects, and 
we have already performed investigations in this direc- 
tion (yet without stress control) in Refs. [27l [39M4l| . A 
similar strategy was taken by other researchers (see, e.g. 

H). 

Numerical results for both, orientational dynamics and 
rheological properties (i.e., stress-strain curves and re- 
sulting viscosities) are presented for stress-controlled sys- 
tems at a fixed value of the "tumbling" parameter [!,[23| 
(that is, the coupling strength between flow and orienta- 
tional motion) , which seems particularly promising based 
on earlier studies under fixed-strain conditions |l2[ . The 
most interesting results are found at temperatures below 
the isotropic-to-nematic (I-N) phase transition, where 
the control method turns out to be crucial for the dy- 
namics of the system. This concerns both the type of bi- 
furcations and their characteristics. Moreover, we show 
how our control scheme can be used to change chaotic 
states seen in systems without stress control into states 
with regular dynamics, or into stationary states. This 
puts our study in the more general context of control in 
non-linear dynamic systems, where the stabilization of 
steady or periodic fix points by feedback control of suit- 
able dynamic variables is a major issue [43j]. 

The paper is organized as follows. In the subsequent 
sections III Atill Bl we summarize the main ingredients of 
the mesoscopic theory for rod-like particles under shear 
flow. Our control method is introduced in Sec. Ill CI Sec- 
tion |lTT] describes our numerical results, with an empha- 
sis on sheared nematic systems. Finally, in Sec. IIVI we 
summarize our main findings and suggest some future 
directions of research. 



II. THEORY 
A. Order parameter and dynamic equations 

We employ a mesoscopic description of the system, 
where the relevant dynamic variable is the orientational 
order parameter averaged over some volume in space. 
In a sheared liquid crystal, this order parameter corre- 



sponds to the time-dependent, 2nd-rank alignment ten- 
sor a = ^/l5/2('uu), where u describes the molecular 
axis and tt7 indicates the symmetric traceless part of a 
tensor. In the isotropic equilibrium state, all compo- 
nents of a are zero. For uniaxial nematic orientational 
order, a = a-\/3/2('nri) with a = \/5S being a scalar pro- 
portional to the Maier-Saupe order parameter 5, and n 
being the corresponding director. In a sheared system, 
however, the assumption of uniaxial order (which leads to 
the so-called Ericksen-Leslie theory of orientational dy- 
namics [iil . [45l | ) remains appropriate only for small shear 
rates and temperatures deep within the nematic phase. 
For intermediate and large shear rates and, in particu- 
lar, in the vicinity of the I-N phase transition, biaxial 
ordering occurs, and a complete description of the non- 
equilibrium dynamics requires the full tensor, &(t). 

Taking a as the relevant variable, the system's equilib- 
rium behavior (in the absence of shear) is governed by 
the Landau-de Gennes free energy [46[ 

$ = l -A{T) & : a - ^-B (a • a) : a + \c (a : a) 2 , (2.1) 

where the notation " :" stands for the trace over the prod- 
uct of the two symmetric tensors, and "•" indicates con- 
ventional matrix multiplication. In Eq. (|2.1[) . T is the 
temperature which is the driving parameter for an I-N 
transition in a thcrmotropic liquid crystal (for lyotropic 
systems, one would replace T by the concentration of the 
particles). 

Switching on an external shear flow characterized by 
a velocity field v, the alignment tensor becomes a time- 
dependent quantity. Here we employ the approach first 
derived from concepts of irreversible thermodynamics 
[22| and later from a generalized Fokker-Planck equa- 
tion for the orientational distribution function |23l - l25| 
(for similar though not identical approaches, see, e.g., 
[HI S3 ) ■ For homogeneous systems the equation of mo- 
tion reads 



(9a 

dt 



2u> x a-2KT-a + r a 7 1 *(a) 



-v^ — r, (2.2) 



where u> = (1/2 )V x v is the vorticity and T = Vv 
is the strain rate tensor. Specializing to the case of a 
plane Couette flow characterized by a flow in x-direction, 
v = 72/e x , the strain rate tensor simplifies to T = 7'e x e y ' 
and u) = — (l/2)7e z . In these formulae, e Q (a = x,y,z) 
denote unit vectors along the coordinate axes. Further, 
the parameters T ap and r a are relaxation time coefficients, 
and k is a dimcnsionless coefficient. Finally, the (tenso- 
rial) quantity $(a) appearing in Eq. (|2.2[) corresponds 
to the derivative of the free energy (|2.1j) with respect to 
the (non-conserved) order parameter, i.e., <fr = d$/d&. 
In the absence of shear (i.e., for v = u> = T = 0) this 
terms governs the relaxation of a(t) towards its equilib- 
rium value determined by $ = 0. We note that the 
Ericksen-Leslie theory of orientational dynamics Q EH 
follows from the present approach, when the alignment 
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tensor is uniaxial and the effect of the shear flow on the 
magnitude of the order parameter can be disregarded. 
Furthermore, as noted already in the Introduction, the 
present approach can be generalized and has been al- 
ready applied to inhomogeneous alignment states [27], HH- 
l4lj | , central steps being the inclusion of gradient elasticity 
terms in the free energy (see also [42| ) and an appropriate 
description of feedback effects of the orientational motion 
on the flow profile [27l [39l - l4l| . Not surprisingly, however, 
these generalizations of the theory imply a substantial in- 
crease of the numerical effort, which is the reason why we 
stick to homogeneous states in the present work. 

We next consider the rheological properties of the 
sheared liquid crystal, which are characterized by its 
pressure tensor, p [36[. The latter consists of an isotropic 
part involving the hydrostatic pressure p, an antisym- 
metric part, and an symmetric traceless part p which we 
refer to as friction pressure tensor J22|. The negative of 
this latter contribution is the so-called shear stress ten- 
sor, <r = —p. This quantity is of special interest since it 
is directly linked with the orientational dynamics in the 
non-equilibrium system. Specifically, one has 



-2r] iso T + p a i, 



(2.3) 



where 77i SO is the so-called second Newtonian viscosity, 
and H 



Pal = 



P^L ( V2^*(a) - 2«^(~ 



(2.4) 



In Eq. (|2.4j) . p is the density of the rod-like particles, m 
is their mass, and /cb is Boltzmann's constant. In the 
absence of shear, one has <fr(a) = and consequently 
PaT = 0. The friction pressure tensor then reduces to the 
first term in Eq. (|2.3[) . which is present also in (sheared) 
systems of spherical particles. A concept to control the 
shear stress externally is formulated in Sec. Ill CI 



B. Explicit equations of motion 

In the practical analysis it is convenient to use scaled 
variables. Details can be found, e.g., in Refs. (22l.l28l.l29j 
and Hl| . The equilibrium behavior of the system is deter- 
mined by the effective temperature, 6 = 9AC/2B 2 . It fol- 
lows that the (first-order) I-N transition occurs at = 1 
(which corresponds to the I-N coexistence point), and 
that the nematic phase (isotropic) is (meta-)stable for 
temperatures 6 < 9/8 (9 > 0). For the description of the 
non-equilibrium system, times and shear rates are made 
dimensionless with a convenient reference time. The lat- 
ter is chosen as the relaxation time of the alignment at the 
isotropic-nematic coexistence, i.e., r re f = r a (9C/2B 2 ). 
Further, the strength of the coupling between the flow 
and the alignment tensor is characterized by the param- 
eter Ak = — (2/3)v / 3'?"ap/' r a- Microscopically, Ak is re- 
lated to the shape of the particles, which can be char- 
acterized by the axis ratio q (for ellipsoidal particles). 



Specifically, one has H A K = 2/( v / 5a K )(g 2 - l)/(g 2 + 1) 
where o,k = 2B/3C. It follows that Ak = for spherical 
particles (q — 1), whereas Ak > for elongated particles 
(q > 1), which is the case considered here. Finally, the 
scaled analogue of the shear rate 7 is denoted as T = r re f7 
(with T rc f defined above). 

Besides scaling, the other important step towards prac- 
tical solution of the problem is to rewrite the tenso- 
rial equation (|2.2|) for the dynamics of a into a set of 
scalar equations for its five independent components, 
tio, cti, • • • , a>i- This is achieved via an expansion of a 
in a standard tensorial basis [III, |4j|. The resulting 
components are linked by simple expressions to their 
cartesian counterparts a a p (with a,0 = x,y,z), that is, 
a cx — (dxx + Oyy), ai oc a xx - a yy , a 2 oc a xy , a 3 cx a xz , 
and ci4 cx a yz [l3|, |49(. Thus, ao (with negative values), 
ai, and ai describe ordering within the shear (x-y) plane, 
whereas 03 and CI4 describe out-of-planc ordering. The 
ordinary (first-order) differential equations for the new 
components of a are given by 



a = 


-00 




Ta 2 


Oi = 


-01 


+ Ta 2 




»2 = 


-02 


-Tai + 


^a r 
— A K r 


a-3 = 


-03 




f 1)04 


(24 = 


-04 




- 1)03, 



"0 



(2.5) 



where the quantities 0o, . . . , 04, which correspond to the 
(dimensionless) components of the free energy derivative 
<fr(a), are given as 



00 = ( 


6 — 'Sag - 


h2a 2 ; 




01 = ( 


+ 6a o - 


h2a 2 ; 


)Ol 


02 = ( 


9 + 6a ~ 


h2a 2 ; 


)a 2 


03 = ( 


9 — 3ao - 


h2a 2 ; 


)a 3 


04 = ( 


9 — 3ao - 


h2a 2 ; 


1 04 



In Eqs. ()2.6|) we have introduced the abbreviation a 2 = 

As outlined before in Sec. Ill AI the rheological proper- 
ties of the system are described by the friction tensor or, 
equivalently, the stress tensor, er. Introducing the ref- 
erence shear modulus G a i, we find from Eqs. (|2.3[) and 

man m 

Gai a = -P = 2 77 iso r -Pai = 277isor + V2G al S al , (2.7) 
where 



^-V6a- # 

oAk 



(2. 
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Here we are particularly interested in the shear stress, 
that is, the x-y component of <r. Expanding a and S al 
into basis tensors one obtains from Eq. (|2.7j) [36j 

o- xy = ?7isor + £ 2 , (2.9) 

where 

2 

S 2 = —7=7— [02 - K(a 2 ^>o + a o 2 )] 
voAk 

-K(a 4 3 + a 3 4 ) (2-10) 
with k = 2k/(3Ak)- 



(at controlled stress conditions) or er xy (controlled shear 
rate), respectively, over an appropriate number of time 
steps. To this end we disregard the initial transient be- 
havior, which may strongly depend on the start condi- 
tions. The overall number of time steps required to ob- 
tain reliable averages strongly depends on the system pa- 
rameters considered. In the present work we consider 
a variety of dynamic "states" characterized either by 
steady-state solutions for a (see Sees. IIII Al and llTl Bp and 
by time-dependent orientational dynamics (Sec. IIII C[) . 
The corresponding total run times vary between 500 and 
150000 time steps. 



C. Controlling the shear stress 

Previous applications of the dynamical equations (|2.5p 
have focussed on the system's orientational dynamics and 
rhcological behavior at constant shear rate, T. The main 
goal of the present work is to explore the "reverse" sit- 
uation where, instead of T, the shear stress is held fixed 
at an (externally imposed) value cr™ p - Indeed, Theologi- 
cal experiments of complex fluids often involve the shear 
stress as an external control parameter rather than the 
shear rate [33l - l35| . 

Within our theory, the instantaneous value of <7 X y is 
determined by V and the instantaneous values for the 
alignment components ak(t) [see Eqs. (|2.9|) and (12.101) ]. 
In order to control the shear stress we follow the exper- 
imental strategy and adjust the shear rate, that is, T 
becomes a time-dependent quantity. This adjustment is 
achieved such that we supplement the five differential 
equations for ak(t) by an additional equation for T(t). 
Specifically, 

Ut( '/ISO 

Here cr xy is the instantaneous shear stress as given by 
(|2.9p . n iso is the (second Newtonian) viscosity discussed 
above, and r g is a control (or "delay") time determining 
the speed of shear stress control. 

D. Numerical solution 

In order to solve numerically the six-dimensional dy- 
namical system consisting of Eqs. (|2.5|) for a-k(t) (k — 
0, . . . , 4) and Eq. (|2.11[) for the shear rate T(t), we employ 
a standard four-step Rungc-Kutta algorithm with adap- 
tive step size control. The resulting time steps are in the 
range At = 0.001 - 0.5. The initial values for ao, ■ ■ ■ , 04 
are set to small, yet unzero values (between 0.01 and 0.1) 
corresponding to a state with weak orientational order. 
The initial shear rate is set to T(0) = 0.01, if not stated 
otherwise. 

In analogy to rheological experiments, which are typ- 
ically not time-resolved, we average the results for T 



III. DYNAMICS AT CONTROLLED SHEAR 
STRESS 

In the following paragraphs we present numerical re- 
sults obtained at various dimcnsionlcss temperatures 
8 around the I-N transition of the equilibrium (un- 
sheared) system. The coupling parameter Ak appear- 
ing in Eqs. (|2.5[) is set to 1.25. This choice is motivated 
by previous investigations at constant shear rate [HI, [l3j 
where, at Ak = 1.25, a variety of dynamics "states" in- 
cluding stationary, oscillatory, and chaotic solutions have 
been found. Here we include such constant-r results as a 
reference. For the constant-stress calculations we choose 
Visa = 0.1 and r g = 1 [see Eq. (|2.11|) ]. if not stated oth- 
erwise. Finally, the parameter k in Eqs. (|2.5[) is set to 
zero. Indeed, earlier studies at constant T indicate that 
the actual value of k has only minor importance for the 
dynamic behavior |36| . 



A. Isotropic phase 

We start by considering the temperature 8 = 1.75 
which is well within the isotropic phase of the unsheared 
system (i.e., a|r=o = 0). Under such conditions, applica- 
tion of an external shear flow yields a stationary, uniaxial 
alignment of the liquid-crystal particles in the x-y plane 
(that is, the flow plane) of the system. Within our the- 
oretical description, shear-induced alignment (indicated 
by an "A" in the following) is reflected by the quanti- 
ties ao, ai and a 2 approaching non-zero, constant values, 
whereas 03 = 04 = 0. The corresponding "constitutive" 
curve <7 xy (r) is shown in Fig.QJ,), where we have included 
data sets from both, constant strain- and constant-stress 
calculations. In both cases, the system quickly reaches 
(after a few hundred time steps) a unique, stationary 
state, independent of the initial conditions. Moreover, 
the results from the two types of calculations are quan- 
titatively consistent to a very high accuracy. Specifi- 
cally, one observes from Fig. []Ji) a monotonic, contin- 
uous, and non-linear (i.e., non- Newtonian) increase of 
the stress with increasing shear rate (note the double- 
logarithmic scale). The corresponding (dimcnsionlcss) 
viscosity i] = d xy /r is plotted in Fig.QJi). The decrease of 
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FIG. 1: a) Shear stress cr xy as a function of the shear rate V 
("constitutive curve") at temperature 9 = 1.75, correspond- 
ing to the isotropic phase of the equilibrium system. Results 
from calculations at fixed Y (cr xy ) are indicated by open circles 
(filled triangles). Parts b) and c) show the corresponding be- 
havior of the viscosity rj and the Maier-Saupe order parameter 
S, respectively. 



7] from its equilibrium value (?yr=o=l-l) to its high-shear 
limit (^r->oo = Vise, = 0.1) reflects a strongly pronounced 
shear-thinning, which is quite typical for complex fluids 
sheared within the isotropic phase Q . Finally, the degree 
of shear-induced (uniaxial) ordering as measured by the 
Maier-Saupe parameter S(T) is plotted in Fig.[T]:). With 
increasing S the corresponding director aligns more and 
more along the y-direction, that is, the direction deter- 
mined by the flow velocity. 

From a more technical point of view it is interesting 
to note that, at the temperature considered, the results 
obtained at fixed <r xy are practically independent of the 
choice of the control time chosen in Eq. (|2.11[) (as tested 
for values of r g between 0.1 and 100). The precise value of 
r g merely affects the amount of time required to reach the 
stationary limit. This is illustrated in Fig. [5^), where we 
plot the time-dependence of the instantaneous shear rate 
F(i) at fixed, imposed stress ct™ p = 0-5 f° r various delay 
times r g . The larger r g , the slower is the approach of T(t) 
towards its asymptotic value ((r) 0.79), as one might 
have expected from the structure of the control equa- 
tion (|2.11|) . The corresponding time dependence of the 
instantaneous shear stress cr xy (i) is plotted in Fig. (2b)- It 
is seen that cr xy (i) relaxes somewhat faster as compared 
to T(t), independent of the actual value of T g . 



B. Shear- induced isotropic-nematic transition 

We next consider the temperature 9 = 1.25. The cor- 
responding equilibrium system is still globally isotropic, 
but much closer to the first-order I-N transition (occur- 
ring at 9 = 1) than the state considered in Scc. lIII Al As a 
consequence of the nearby equilibrium transition, shear- 
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FIG. 2: Instantaneous shear rate (a) and shear stress (b) as 
functions of time under controlled stress conditions (er™ 13 = 
0.5) at 9 — 1.75. Included are results for four relaxation times 
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FIG. 3: Same as Fig. [TJ but at temperature 9 = 1.25 close 
above the I-N transition of the equilibrium system. 



ing induces a non- equilibrium transition from a "parane- 
matic" state characterized by weak orientational ordering 
into a state with pronounced stationary and (almost) uni- 
axial nematic order. We note that the equilibrium system 
has a mctastable nematic phase only for lower tempera- 
tures, that is, for 9 < 9/8 = 1.125). The shear-induced 
transition is illustrated in Figs.[3Ji)-c), where we plot the 
stress-strain relation at 9 = 1.25 together with corre- 
sponding results for the viscosity and the Maier-Saupe 
parameter. For imposed shear rate, the non-equilibrium 
transition occurs T w 0.178 where we observe pronounced 
discontinuities in all of the three quantities investigated. 
Specifically, the jump-like increase of S is accompanied 
by jumps in <7 xy (and, consequently, 77) towards smaller 
values, indicating a sudden shear-thinning of the system. 
We note that the occurrence of the shear-induced tran- 
sition has been realized, on a theoretical basis, already 
quite some time ago jlOl lll| . A recent experimental re- 
alization involves colloidal suspensions of tobacco viruses 
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FIG. 4: Shear stress cr xy as a function of the shear rate F at 
temperature 8 = 0, within the nematic phase of the equilib- 
rium system. Results from calculations at fixed V (a xy ) are 
indicated by open (filled) circles. Calculations at fixed stress 
have been performed at r g = 1. Fhe abbreviations "KF", 
"W", "C", and "A" indicate the orientational behavior and 
are explained in the main text. Fhe inset shows the corre- 
sponding behavior of the viscosity. 
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Coming back to our model system, the data in Fig. [3Ja,) 
imply that there arc values of er xy where two dynamic 
states characterized by different shear rates T can " coex- 
ist" . Interestingly, when controlling the shear stress [see 
black triangles in Fig. Ek)], the system remains in the 
branch with lower T up to largest value of cr xy where co- 
existence is possible. This observations turns out to be 
essentially independent of the details of the stress con- 
trol, i.e. the control time r g (and the initial value for T). 
This rather robust behavior is in sharp contrast to the 
non-equilibrium transitions at lower temperatures to be 
discussed next. 



C. Shearing within the nematic phase 

As a representative example for the system's low- 
temperature dynamics we consider the temperature 9 — 
0. Here, the nematic state is globally stable already 
at zero shear, and we are at the lower limit of (meta- 
)stability of the isotropic phase. Our results for the 
shear stress as function of T are displayed in Fig. [4] 
All calculations have been started from almost random 
initial states (see Sec. Ill D[) . Earlier numerical studies 
fl2j with controlled shear rate and Ak = F25 have al- 
ready revealed the occurrence of various types of time- 
dependent dynamics of the nematic director. These in- 
clude " symmetry-breaking" states where the director ro- 
tates out of the shear plane (i.e., a (t) > 0, (t) ^ 0, 
<24 (i) 7^ 0). Our present results at constant T are fully 
consistent with those in [l2j • Specifically, we observe for 
low shear rates (r < 1.7) so-called kayaking-tumbling 
(KT) motion, where director performs out-of-plane os- 
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FIG. 5: Enlarged version of the stress-strain curve in Fig. [4] 
The lines indicate the result of constant-stress calculations 
with different values of the relaxation time r E . 



cillations, and its projection onto the shear plane de- 
scribes an ellipse (i.e., tumbling). Corresponding orbits 
illustrating the time dependence of 01-04 are shown in 
Figs. [5] and H] Increasing T there exists a (small) win- 
dow (1.7 < r < 1.9) characterized by kayaking-wagging 
(KW) motion, where the director performs out-of-plane 
oscillations while its projection in the shear plane de- 
scribes wagging. From Fig. 2] we see that the transition 
into the KW state is supplemented by a sudden decrease 
of the shear stress. Upon further increase of T we are 
re-enterting the KT state until T w 3.5. In the subse- 
quent range of shear rates the system has chaotic (C) fea- 
tures, with the long-time dynamics being strongly depen- 
dent on the initial conditions and the largest Fyaponov 
exponent being positive (fo r detailed discussion of the 
chaotic regime, see Refs. |I3.ll3j ). Finally, for shear rates 
r > 4.15 there is a further transition into a flow-aligned 
(A) state where the director is arrested to a constant 
direction within the shear plane. 

Remarkably, not all of these states are observable when 
we take the shear stress as a control variable. This can 
be seen from the data points indicated by (black) tri- 
angles in Fig. [4] Specifically, our controllcd-stress cal- 
culations do neither yield the wagging motion, nor the 
complex, chaotic behavior. Instead, the system jumps 
directly from the KT into the shear-aligned (A) state. 
Moreover, it turns out that the position of this jump in 
the r-er xy plane can be manipulated, to some extent, by 
the control time r g (indeed, the initial value of the shear 
rate has much less effect). The role of r g is illustrated 
in Fig. [5l where we present results for various values of 
T g on an enlarged scale. The larger the control time, 
the longer the system remains on the low-r branch. In 
this sense, an increase of r g stabilizes the periodic KT- 
state. Interestingly, the chaotic region (C) does not ap- 
pear in stress-controlled calculations even for the largest 
values of r g considered. Indeed, as we will demonstrate 
in Sec. IIIID1 a sudden onset of shear-stress control can 
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FIG. 8: Phase portraits of the in-plane components of the 
pressure tensor, ei2 versus a\, for different control times r g . 
The average shear stress is fixed at cr™ p = 0.575. 



FIG. 6: Instantaneous stress as function of time for different 
control times r g . The average shear stress is fixed at er™ p — 
0.575. The resnltinp' time stens are At. 0.01 for all values 
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FIG. 9: Same as Fig. [8l but for the out-of plane components 
d4 versus a^. 



FIG. 7: Instantaneous shear rate as function of time for dif- 
ferent control times r g . Parameters as in Fig. HJ] 



actually stabilize these chaotic states. 

Even away from the KT-A transition the dynamics 
in the stress-controlled, ncmatic system depends sig- 
nificantly on the control time, r g . This is illustrated 
in the subsequent figures IBJIHl where we consider vari- 
ous time-dependent quantities at an imposed stress of 
ct™ p = 0.575. According to Fig. 21 this value corre- 
sponds to a state within the KT regime (and below the 
KT-A transition). To start with, the time-dependence 
of the instantaneous stress o~ xy (t) and that of the shear 
rate T(t) are displayed in Figs. |6] and [3 respectively. 
For all control times considered, both quantities display 
pronounced oscillations even in the long-time limit. In 
particular, o~ xy (t) reaches its externally imposed value of 
°xy P = 0.575 only on the average [see Fig. [6]. The ampli- 
tude of the oscillations strongly depends on r g . Specifi- 
cally, an increase of the control times tends to damp out 
the oscillations. Similar behavior is seen in Fig. [7]for the 
shear rate 

The oricntational dynamics corresponding to the KT 



state at cr™ p = 0.575 is illustrated in Figs. [5] and H 
where we show phase portraits (orbits) of both, the in- 
plane motion (02 as function of ai) and the out-of-planc 
motion (04 as function of 03). The ellipse-shaped limit 
cycle characterizing the functions 02(01) [see Fig. [5] is 
typical for the KT motion, where the nematic director 
displays tumbling motion in the shear plane. A closer 
analysis of the time-dependence of the individual com- 
ponents reveals that the transient behavior observed in 
the period between rj 400 and w 1400 time steps (i.e., 
well before the asymptotic KT state is reached) some- 
what depends on the control time, r g . Specifically, at 
T g = 1 there is transient flow alignment where ao, &i, 
a2 are constant, and 03 ~ 0, 04 w 0. At r g = 10 one 
finds instead a short period of in-plane tumbling, where 
Oi and 02 become time-dependent. These differences in 
the alignment dynamics are mirrored by corresponding 
differences in the transient behavior of the instantaneous 
shear rate and stress, respectively (see insets in Figs. [5] 
and [7]) . At larger times the orbits related to different 
values of the control time r g look rather similar. This is 
in accordance to earlier findings on the robustness of the 
orientational dynamics with respect to time-dependent 
perturbations of the shear rate [501 ] - 
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FIG. 10: Shear stress and shear rate as functions of time in 
systems where the control method switches from fixed shear 
rate (T = 3.7) to fixed shear stress (ct™ p = 2.6, r g = 1). The 
switch occurs at a) N s — 5 x 10 4 time steps, b) N B = 1 x 10° 
time steps. The resulting time steps are At « 0.01 for both 
values of N s . 



D. Stabilization of chaotic states 

An interesting feature of our results within the nematic 
phase is that under controlled stress conditions, the tran- 
sition from the KT into the shear-aligned state occurs 
before the system enters the chaotic (C) regime observed 
in calculations at fixed V (sec Fig. [4j) . This suggests that 
stress control might have a stabilizing influence on shear- 
rate-induced chaotic states. To explore this conjecture 
we have performed test "experiments" where the con- 
trol method was switched from fixed shear rate to fixed 
stress after a predefined number N s of time steps. That 
is, for N < N s the oricntational dynamics was calcu- 
lated via the five equations (|2.5|) . whereas at N > N s 
these equations were used together with Eq. (|2.11|) . As 
an representative example, we consider the behavior ob- 
tained at T = 3.7 (N < N s ) and a™ p = 2.6 (N > N s ). 
Numerical results for the instantaneous shear stress and 
-rate, are displayed in Fig. 1101 where we have included 
data for two values of 7Y 8 . The corresponding behavior 
of the in-plane components of the director, a\ and 02, is 
shown in Fig. 111! In both figures, one clearly recognizes 
a drastic change of the dynamics upon the switch from 
shear- to stress control. Before the switch one observes 
highly irregular dynamics of all non-fixed quantities, that 
is cr xy (<), ai(t), a,2{t) (and the same holds for the other 
components of the tensor a). Closer results reveals that 
the system at fixed V = 3.7 is indeed in a chaotic state 
[HI, [Tj]. After the switch, that is, under controlled stress, 
the dynamics becomes regular. Interestingly, the charac- 
ter of this final state depends on both, the starting point 
(N B ) of the stress control, and the corresponding delay 
time r g . For the choices N s = 5 X 10 4 and r g = 1 [see 
Figs. fTDk) andlTTa)]. the systems displays KT dynamics 
characterized by oscillatory (yet regular) motion of all 



FIG. 11: Same as Fig. 1101 but for the components ai and 02 
of the alignment tensor. 



components a\-ai and corresponding oscillations of T(t). 
The average shear rate is (r) w 3.1, which is again a 
typical value for the KT regime as seen in Fig. |4j Per- 
forming then the switch at the later time N s = 1 x 10 5 
(and T g = 1) [see FigsfTUb) andHTb)]. all of the displayed 
quantities become constant. This signals that we are en- 
tering a shear aligned (A) state, characterized by a shear 
rate of (r) rj 25.8. It is noteworthy that the same A 
state is achievable with the choices N s = 5 x 10 4 and 
r g = 10 (or 100). Taken altogether, these results reveal 
that stress control can indeed stabilize the dynamics in 
systems characterized by unstable motion without con- 
trol. Moreover, the character of the final state is tunable 
by the details of the control scheme, that is in our case, 
the starting time of the stress control, and the corre- 
sponding delay time. 



IV. CONCLUSIONS 

In this article we have proposed a mesoscopic theory 
to investigate the oricntational dynamics and rheology 
of sheared complex fluids at fixed, externally imposed 
shear stress cr xy . Our strategy to control er xy consists 
of supplementing the previously established equations of 
motion for the orientational order parameters [l3T ] by a 
relaxation equation for the shear rate, T. The speed of 
this relaxation can be tuned by a control ("delay") time. 
Based on the resulting six-dimensional dynamical sys- 
tem we have numerically investigated sheared fluids of 
rod-like particles at several temperatures around the I- 
N transition at T = 0. At high temperatures, where 
the equilibrium system is isotropic, shear stress control 
yields the same asymptotic behavior (that is, flow align- 
ment) as that encountered at fixed T. Within the low- 
temperature nematic phase, however, the control vari- 
able (stress or strain) influences both, the location and 
the type of bifurcations between different dynamic states. 
This finding is relevant not only from a theoretical point 
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of view but also for experiments; indeed, in rheological 
set-ups it is often the shear stress rather than the shear 
rate which can be fixed externally. In our model, the 
stress-controlled nematic fluid is characterized by only 
one non-equilibrium transition from a periodic (kayaking- 
tumbling) into a flow-aligned state, contrary to its T- 
controllcd counterpart which displays, in addition, wag- 
ging and chaotic motion. Furthermore, the dynamics at 
fixed stress is strongly influenced by the delay time r g . 
Large values of the latter tend to damp oscillatory mo- 
tion. Moreover, the delay time can be used as a parame- 
ter to tune the position of the system's transition from a 
low-r to a high-r branch. Finally, we have demonstrated 
that our control scheme is capable of stabilizing regular 
(periodic or stationary) states out of chaotic orbits ap- 
pearing at fixed T. 

In the present paper we have restricted ourselves to 
the investigation of spatially homogeneous states. This 
is clearly an oversimplification in view of experiments 
pointing to shear-induced inhomogeneities (see, e.g., 
[2ll|37|,|38|). However, as mentioned before, our approach 
is in principle fully capable of treating inhomogeneities 
as demonstrated, e.g., in Refs. 27 1 and [HI], the latter 
study dealing with a full three-dimensional spatially re- 
solved flow problem. Interestingly, there is a recent study 
of Das et al. J2(| suggesting that the chaotic orientational 
behavior of our present system at constant T transforms 
into a spatiotemporal chaos, when spatial fluctuations are 
allowed in the theory (via gradient terms in the free en- 
ergy). This suggests that the stress-controlled dynamics 
reported in present work may also be spatially extended. 
We note, however, that a full description of the spatially 
resolved problem would include feedback effects of the ori- 
entational motion on the flow [27J (which are neglected in 
[2r| . where V was assumed to be a constant). Given the 
complexity of the resulting numerical study we consider 
the present findings as first indicators for the true dy- 
namics and as useful starting points for a more detailed 
investigation. A further interesting direction would be 



an extension of the present study towards polar systems 

Finally, the findings in the present study, particularly 
those regarding chaotic states, are also interesting in the 
broader context of control in non-linear dynamic systems 
[43I ]. In particular, there is currently an intense research 
on systems such as lasers [52j], neural systems [53|, ex- 
citable media [HI, [55|, and Brownian ratchet systems 
[5rj [57j , where unstable dynamic behavior is controlled 
via time- delayed feedback control (TDAS) [EH]. In this 
framework the equations of motion are supplemented by 
" control force" terms involving differences of an appropri- 
ate system variable (or output quantity) at a time t and 
an earlier time t—r. The control force thus vanishes when 
a specified stationary or periodic state is reached. This 
strategy is somewhat different from the present approach 
where the control is realized by an additional equation of 
motion for the shear rate T, which is varied until cr xy 
approaches an externally imposed value. Despite these 
differences, a common feature of our approach and the 
TDAS method is that the time parameter (r or r g , re- 
spectively) plays a crucial role for the asymptotic dynam- 
ics (in fact, the present scheme could be viewed as feed- 
back control with additional low-pass filtering involving 
an exponentially distributed delay [H^]). Indeed, in our 
model the actual value of T g decides which type of non- 
chaotic state the unstable system selects when starting 
the control within the chaotic regime. These similarities 
suggest that control concepts established in non-linear 
systems, such as TDAS, may be very fruitful not only 
for a deeper understanding of the rheological properties, 
but also for the orientational dynamics. Work in this 
direction is under way. 
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